Modeling the quantum interference signatures of the Ba II D2 
4554 A line in the second solar spectrum 



H. N. Smitha 1 , K. N. Nagendra 1 , J. 0. Stenflo 2 ' 3 and M. Sampoorna 1 

1 Indian Institute of Astrophysics, Koramangala, Bangalore, India 
2 Institute of Astronomy, ETH Zurich, CH-8093 Zurich, Switzerland 
3 Istituto Ricerche Solari Locarno, Via Patocchi, 6605 Locarno-Monti, Switzerland 
smithahnOiiap . res . in; knn@iiap . res . in; stenflo@astro.phys.ethz.ch; sampoornaOiiap . res . in 



ABSTRACT 

Quantum interference effects play a vital role in shaping the linear polarization profiles of solar 
spectral lines. The Ba II D2 line at 4554 A is a prominent example, where the F-state interference 
effects due to the odd isotopes produce polarization profiles, which are very different from those 
of the even isotopes that have no F-state interference. It is therefore necessary to account for the 
contributions from the different isotopes to understand the observed linear polarization profiles 
of this line. Here we do radiative transfer modeling with partial frequency redistribution (PRD) 
of such observations while accounting for the interference effects and isotope composition. The 
Ba II D2 polarization profile is found to be strongly governed by the PRD mechanism. We show 
how a full PRD treatment succeeds in reproducing the observations, while complete frequency 
redistribution (CRD) alone fails to produce polarization profiles that have any resemblance with 
the observed ones. However, we also find that the line center polarization is sensitive to the 
temperature structure of the model atmosphere. To obtain a good fit to the line center peak of 
the observed Stokes Q/I profile, a small modification of the FALX model atmosphere is needed, 
by lowering the temperature in the line-forming layers. Because of the pronounced temperature 
sensitivity of the Ba II D2 line it may not be a suitable tool for Hanle magnetic-field diagnostics 
of the solar chromosphere, because there is currently no straightforward way to separate the 
temperature and magnetic-field effects from each other. 
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The linearly polarized spectrum of the Sun 
known as the Second Solar Spectrum exhibits sig- 
natures of a number of physical processes not seen 
in the intensity spectrum and which may also be 
used to diagnose weak and turbulent magnetic 
fields that are inaccessible with the ordinary Zee- 
man effect. Many of the most prominent spectral 
lines in the Second Solar Spectrum, like Nai Di 
and D2, Ball Di and D2, and Call H and K, 
are governed by effects of quantum interference 
between states of different total angular momen- 
tum (J or F states). The profound importance 



of quantum interference for the formation of these 
lines was first demo nstrated both observa tionally 
and t heoretically bv lStenfio ( 1980l see also lStenflol 
1997I K 

Quantum interference between the J-states 
governs the linear polarization signatures of 
atomic multiplets such as the Call H and K, Nai 
Di and D2, and Cr 1 triplet lines. The theory 
of such J-state interference was developed to in- 
clude the e ffects of partia l frequ e ncy re distribution 
(PRD) by ISmitha et all (|201ll . |2013h . This the- 
ory was then applied successfully to model the 
linear polarizati on profiles of the Cr 1 triplet at 
5204-5208 A bv lSmitha etaD (|2012al ). 
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When an atom possesses nuclear spin the J 
states are split into F states (hyperfine structure 
splitting). The quantum interference between the 
F states produces depolarization in the line core. 
Examples of lines governed by F-state interference 
are the Nai D 2 , Ban D 2 , and Sen line at 4247 A. 
The Ba n D 2 line is due to the transition between 
the upper fine structure level J = 3/2 and the 
lower level J = 1/2 (see Figure QJa)). In the odd 
isotopes of Ba, both the upper and lower levels 
undergo hyperfine structure splitting (HFS) due 
to the nuclear spin I s = 3/2, resulting in four up- 
per and two lower F-states (see Figure H{b)). The 
quantum interference between the upper instates 
needs to be taken into account in the modeling 
of the Ball D 2 line. The odd isotopes contribute 
about 18% of the total Ba abundan ce in the Sun 
(c.f. Table 3 of lAsplund et al.ll2009h . The remain- 
ing 82% comes from the even isotopes, which are 
not subject to HFS (because I s = 0). 

The intensity profile of the Ba II D 2 line has 
earlier been studied exten sivel y, for exampl e by 
Holweeer fc Muellerl (|l974T ) and|Rutten| (|l978l and 
the references cited therein). Some of these stud- 
ies aimed at determining the solar abundance of 
Ba. Observations with th e high precision spectro- 
polarimeter ZIMPOL by IStenflo fc Kelierl (jl997l ) 
clearly revealed the existence of three distinct 
peaks in the linear polarization (Q/I) profiles of 
the Ba n D 2 line. The nature of these pea ks could 
subseq uently be theoretically clarified by IStenflo 
(|l997l ). who used the last scattering approxima- 
tion to model the Q/I profiles. It was demon- 
strated that the central Q/I peak is due to the 
even isotopes of Ba, while the two side peaks are 
due to the odd isotopes. 

Using a similar last scattering approximation, 
the mag netic sensitivity of t he Ba line was ex- 
plored by [BeUuzzTetjiL ( 2007 ). Both these papers 
however did not account for radiative transfer or 
PRD effects. The potential of using the Ban D 2 
line as a diagnostic tool for chromospheric weak 
turbulent magnetic fields and the important rol e 
of PRD were discussed bv lFaurobert et al~ (2009), 
but the treatment was limited to the even Ba 
isotopes, for which HFS is absent. In contrast, 
our radiative-transfer treatment with PRD in the 
present paper includes both even and odd isotopes 
and the full effects of HFS with F-state intcrfer- 



The theory of F-state interference with PRD in 
the non-mag netic collisionle s s regim e was recently 



developed in ISmitha et al.l (|2012bl . herafter PI). 



The PRD matrix was also incorporated into the 
polarized radiative transfer equation. The trans- 
fer equation was then solved for the case of con- 
stant property isothermal atmospheric slabs. In 
the present paper we extend the work of PI to 
solve the line formation problem in realistic 1-D 
model atmospheres in order to model the Ba n D 2 
line profile observed in a quiet region close to the 
solar limb. 

The outline of the paper is as follows: In Sec- 
tion [5] we present the polarized radiative transfer 
equation which is suitably modified to handle sev- 
eral isotopes of Ba. In Section [3] we present the 
details of the observations. In Section @] we dis- 
cuss the model atom and the model atmosphere 
used. The results are presented in Section [5] with 
concluding remarks in Section [5] 

2. Polarized line transfer equation with in- 
state interference 

The polarization of the radiation field is in 
general represented by the full Stokes vector 
(I, Q, U, V) T . However, in the absence of a mag- 
netic field Stokes U and V are zero in an axisym- 
metric 1-D atmosphere. Hence in a non-magnetic 
medium the Stokes vector (J, Q) T is sufficient to 
express the polarization state of the radiation field. 
The trans fer equation in the re duced Stokes vector 
basis (see ISmitha et al. 2012a ) is 



dl (\,n, z) 

dz 



-k tot (\,z)[I(\,n,z)-S(\,z)}, (1) 



with positive Q defined to represent linear polar- 
ization oriented parallel to the solar limb. The 
quantities appearing in Equation (p} arc defined in 
the reference mentioned above. However, we need 
to generalize the previous definitions of opacity 
and source vector to handle even and odd isotope 
contributions together. 

The total opacity k to t(X,z) = ki{z)<f> g {\,z) + 
cr c (A, z) + fcth(A, z), where a c and fc t h are the con- 
tinuum scattering and continuum absorption coef- 
ficients, respectively. In the present treatment, fcth 
also includes the contribution from the blend lines, 
which are assumed to be depolarizing and hence 
are treated in LTE. fc; is the wavelength averaged 
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6s 2 S 1/2 ;I s =3/2 



- F=2 



6s 2 S„ 



Fig. 1. — (a) represents the Ba n model atom for the even isotopes, while for the odd isotopes the atomic 
model is modified by replacing two of the levels, 2 P^/2 an d 2 <Si/2j with their hyperfine structure components, 
as shown in (b). In (b), I s is the nuclear spin. The energy levels are not drawn to scale. 



absorption coefficient for the J a — > J& transition. 
J a and Jb are the electronic angular momentum 
quantum numbers of the lower and upper level, 
respectively. <j) g is the Voigt profile function writ- 
ten as 

<j) g (\, z) = 0.822 cf> e (X, z) + 0.178 O (A, z). (2) 

4> e (X, z) is the Voigt profile function for the even 
isotopes of Ba n corresponding to the J a = 1/2 — » 
Jb = 3/2 transition in the absence of HFS. The 
profile function for the odd isotopes is <j> (X,z), 
which is the weighted sum of the individual Voigt 
profiles (j)(\F b F a ,z) representing each of the F a —> 
Fb absorption transitions. Here F a and Fb are the 
total angular momentum quantum numbers of the 
initial and the intermediate hyperfine split levels, 
respectively. <fi (X,z) is the same as 0hfs(A,z) 
defined in Equation (7) of PI and is given by 



2 5 
—<f>(\oi,z) + —^(Xn,z) 



5 1 

+— 4>(Xii,z) + — <j>{X 12 ,z) 

5 14 

+ ^{M2,Z) + —<t>{X32,z) 



(3) 



The 17.8% of 4> (X, z) in Equation @ contains 
contributions from both the 135 Ba (6.6%) and 
137 Ba (11.2%) odd isotopes. 

The reduced total source vector S(X, z) appear- 



ing in Equation ([T]) is defined as 

c ,, s _ ki(z)Si(X,z) +a c (X,z)S c (X,z) 

fc th (A, z)5 t h(A, z) + e ki(z)(pg(X, z)<S t h(A, z) 



k to t(X,z) 



(4) 



for a two-level atom with an unpolarized lower 
level. Fo r the c ase of Ba n D2 it was shown by 
Dero uich (2008) that any ground level polariza- 



tion would be destroyed by elastic collisions wi th 
hydrogen atoms (see also iFaurobert et ah 20091) . 

In Equation flD, <Sth = (B\,0) T , where B\ is 
the Planck function. <S/(A, z) is the line source 
vector defined as 



Si(X,z) 



+00 ^ 
2j_ 



^(A,A',z)*( M ') 
xX(X',fx',z)dfi' dX', (5) 



with 



fc(X,X',z) = 0.822 7?. e (A,A',z)+0.178 7?. o (A,A',z) 

(6) 

Here TZ (X,X',z) is a (2 x 2) diagonal matrix, 
which includes the effects of HFS for the odd iso- 
topes. Its elements are H = diag (R®, 71%), where 
TZ^ are the redistribution function components 
for the multipolar index K, cont aining both type - 
II and type-Ill redistribution of iHummerl (|1962n . 
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The expression for 1Z K is obtained by the quantum 
number replaceme nt L — » J; J — > F; S —> I s 



in Equation (7) of ISmitha et al.l (|2012bl . see also 



Smitha et al. 20131 and PI). In our present com- 



putations, we replace the type-Ill redistribution 
functions by CRD functions. We have verified 
that bot h of these give nearly identical re sults 
(see also lMihalaslll978t ISmitha et al]l2012ah and 
such a replacement drastically reduces the com- 
putation time. The redistribution matrix for the 
17.8% of the odd isotopes includes the contribu- 
tions from the individual redistribution matrices 
for the 135 Ba and 137 Ba isotopes. 

1Z e (X, A', z) is also a (2 x 2) diagonal matrix for 
the even isotopes without HFS. Its elements IZ 1 ^ 
are the redistribution functions corresponding to 
the J a = 1/2 -> J b = 3/2 ->• J f = 1/2 scatter- 
ing transition. They are obtained by setting the 
nuclear spin I s = in H (X, X',z). An expression 
for 7t P {X, X' ,z ) can be found inlDomke k, Hubenv 
(1988] ) and in Bommierl ( 19971 see also iNagendra 
Il994l ISampoornal 20 111 ) in the Stokes vector basis. 
It is the angle averaged versions of these quantities 
that are used in our pre sent computat i ons. A s has 
been demonstrated in ISupriva et al.l (|2013j), the 
use of the angle-averaged redistribution matrix is 
sufficiently accurate for all practical purposes. 

Like in ISmitha et all (|2012bl ) we use the two 
branching ratios defined by 



A = 



BW = 



r fi + Tj 



r, 



D {K) 



(7) 



and Ti are the radiative and inelastic colli- 
sional rates, respectiv ely. T e is the elastic c ollisio n 
rate computed from Barklem fc O'Mara (|l998l ). 

D {K) 

are the depolarizing elastic collision rates 
with p(°) = The is compute d using (see 
Derouichll2005 iFaurobert et al.ll2009h 



= 6.82 x 10 - 9 7i H (T/5000) - 40 
+7.44 x 10" 9 (l/2) 1 - 5 n H (T/5000)°- 38 exp(A£;/fcT), 

(8) 

where «h is the neutral hydrogen number density, 
T the temperature, and AE the energy difference 
between the 2 Pi/2 and 2 i"3/2 fine structure lev- 
els. In the present treatment we neglect the col- 
lisional coupling between the 2 P^/2 level and the 



metastable 2 D 5 / 2 level. The importance of such 
collisions for the line cente r polarization of Ba 11 
D2 has been pointed out bv lDerouich (2008). who 
showed that the neglect of such collisions would 
lead to an overestimate of the line core polariza- 
tion by ~ 25%. This in turn would cause the mi- 
croturbulent magnetic field (B tur b) to be overes- 
timate d by ~ 35%, as shown bv iFaurobert et al 
(2009). However, the aim of the present paper 



is not to determine the value of -B tur b but to ex- 
plore the roles of PRD, HFS, quantum interfer- 
ences, and the atmospheric temperature structure 
in the modeling of the triple peak structure of the 
Ba 11 D2 linear polarization profile. 

Frequency coher ent scattering is ass umed in 
the continuum (see ISmitha et al. 2012a ) with its 
source vector given by 



S c (X,z) 



1 



*( fJ ,')Z(X, f i',z)d»'. (9) 



The matrix ^ is the Rayleigh sc attering phas e 
matrix in the reduced basis (see Frisch 2007). 
The line thermalization parameter e is defined by 
e = Ti/(T R + Tj). The Stokes vector (J, Q) T can 
be computed from the irreducible Stokes vector 
X by simple transformations given by (see iFrischl 
20071) 



i(x,e,z) 



Q(x,e,z) 



+^=(3cos 2 6-l)l%(X,n,z), 
(1 - cos 2 6»)Z 2 (A,^,z), (10) 



2^2 



where 9 is the colatitudc of the scattered ray. 
The scatterin g geom etry is shown in Figure 1 of 
Anusha et all (|201lh . 



3. Observational details 

The observed polarization profiles of the Ba II 
D2 line that are used in the present paper for mod- 
eling purposes were acquired by the ETH team 
of Stenflo on June 3, 2008, using their ZIMPOL- 
2 imaging polarimetry system (|Gandorfer et al 



2004) at the THEMIS telescope on Tenerife. Fig- 
ure[2]shows the CCD image of the data recorded at 
the heliographic north pole with the spectrograph 
slit placed parallel to the limb at \i = 0.1. The po- 
larization modulation was done using Ferroelectric 
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Liquid Crystal (FLC) modulators. The spectro- 
graph slit was 1" wide and 70" long on the solar 
disk. The resulting CCD image has 140 pixels in 
the spatial direction and 770 pixels in the spectral 
direction. The effective pixel size was 0.5"spatially 
and 5.93 mA spectrally. The observed profiles 
used to compare with the theoretical ones have 
been obtained by averaging the / and Q/I images 
in Figure [2] over the spatial interval 40"- 52". 

The recording presented in Figure [5] does not 
show much spatial variation along the slit, since it 
represents a very quiet region. However, record- 
ings near magnetic regions made during the same 
observing campaign with ZIMPOL on THEMIS 
exhibit large spatial variations. It has long been 
known that all strong chromospheric scatter- 
ing lines (like the Ca I 4227 A, Na I D 2 , Sr II 
4079 A line, etc) have such spatial variations. Our 
observations confirm that the Ba II D 2 line is no 
exception, which means that it is sensitive to the 
Hanle effect like the other chromospheric lines. 
Observations of spatial variations of this line have 
also been carried out by Lopez Ariste et al. ( 20091 ) 
and lRamelli et all <|2009h . 



4. Modeling procedure 

the polarization profiles of the 



To 
Ba II 



model 
D 2 lin 



the one described in Ho 


zreuter et al. (2005, see 


also Anusha et alJ 2011 


Anusha et al. 2010l 


Smitha et al.l 2012a). It involves the computation 



of the intensity, opacity and collisional rates from 
the PRD-capable MALI (Multi-level A pproximate 
Lamb da Iteration) code developed by lUitenbroek 



(|200lL referred to as the RH-code). The code 
solves the statistical equilibrium equation and 
the unpolarized radiative transfer equation self- 
consistently. The opacities and the collision rates 
thus obtained are kept fixed, while the reduced 
Stokes vector X is computed perturbatively by 
solving the polarized radiative transfer equation 
with the angle-averaged redistribution matrices 
defined in Section [2] 

Such a procedure requires a model atom and a 
model atmosphere as inputs to the RH-code. The 
details of the model atom and the atmosphere are 
discussed in the next subsections. 



4.1. Model atom 

Three different atom models are considered, 
two for the odd and one for the even isotope. The 
atom model for the even isotope ( 138 Ba) is given 
by the five levels of Figure QJa), while for the odd 
isotopes ( 135 Ba and 137 Ba) the model is extended 
to include the hyperfine splitting as described by 
Figure W[b). We neglect the contribution from 
other less abundant even isotopes. The wave- 
lengths of the six hyperfine transitions for the odd 
isotopes are taken from Kurucz' database and are 
listed in Table 1. These transitions are weighted 
with their line strengths given in Equation (|3|) (see 
Table 1). 

4.2. Model atmosphere 

We present the results computed for some of 
the standard realistic 1 -D model atmosphere s, like 
FALA, FALF, FALC (jFontenla et al.l Il9935 and 



FALX ( AvrettJ 1995 ). Among these four models 
FALF is the hottest and FALX the coolest. Their 
temperature structures are shown in the top panel 
of Figure [3J However, as will be discussed below, 
we find that a model atmosphere that is cooler 
than FALX is needed to fit the observed profiles. 
The new model, denoted FALX, is obtained by 
reducing the temperature of the FALX model by 
about 300 K in the height range 500 - 1200 km 
above the photosphere. 

We have verified that such a modification of 
the FALX model does not significantly affect the 
intensity spectra. In contrast, the Q/I spectra 
turn out to be very s ensitive to su c h temp era- 
ture changes. Like in Smitha et al.l ( 2012a ). we 



test the FALX atmosphere by computing the limb 
darkening function for a range of wavelengths and 
fi va lues and compa r e it w ith the observed data 
from iNeckel fc Labs! (1994). This is shown in the 



bottom panel of Figure [3] One can see that FALX 
and the standard FALX fit the observed center-to- 
limb variation equally well. Therefore small modi- 
fications of the temperature structure to achieve a 
good fit to the observed Q/I profile can be made 
without affecting the model constraints imposed 
by the intensity spectrum. 

5. Results 

In the following we discuss the modeling details 
and the need for a model atmosphere that is cooler 
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than FALX. This helps us to evaluate the temper- 
ature sensitivity of the Ba II D2 line and its use- 
fulness for magnetic-field diagnostics. In addition 
we demonstrate the profound role that PRD plays 
for the formation of the polarized line profile. 

5.1. Modeling the Ba 11 D2 line profile 

5.1.1. Modeling details 

From the three Ba II atom models described in 
Section 14.11 we obtain three sets of physical quan- 
tities (two for the odd isotopes and one for the 
even isotope) from the RH code. These quanti- 
ties include line opacity, line emissivity, continuum 
absorption coefficient, continuum emissivity, con- 
tinuum scattering coefficient, and the mean inten- 
sity. The mathematical expressions used to com- 
pute these v arious quant i ties fo r the even isotopes 



are given in lUitenbroekl (|200l[ ) . For the odd iso 



topes, the profile functions in these expressions are 
replaced by </> (A, z) defined by Equation ((3j). 

The three sets of quantities are then combined 
in the ratio of their respective isotope abundances 
and subsequently used as inputs to the polariza- 
tion code. 

5.1.2. Temperature sensitivity 

The polarization profiles thus computed for the 
various model atmospheres are shown in Figure |H 
displayed separately for the even, odd, and com- 
bined even-odd cases in three different panels. The 
Stokes Q/I profiles in Figure 0] are computed by 
setting the total abundance of Ba in the Sun equal 
to the abundance of even isotopes in the first 
panel; the abundance of odd isotopes in the second 
panel; and a fractional abundance of even (82%) 
and odd (18%) isotopes in the third panel. The 
profiles in the first panel can be compared to the 



result s presented in Figure 6 of iFaurobert et al 
(2009). As seen from the first panel, the ampli- 



tude of the central peak for the even isotopes is 
very sensitive to the temperature structure of the 
model atm osphere in contradictio n with the con- 



clusions of Faurobert et al. (2009). Also in their 



paper, the amplitude of the central peak obtained 
from the FALC model is larger than the one ob- 
tained from FALX, which is opposite to our find- 
ings (although it could be that the version of the 
FALC model they used is not identical to the one 
that we have used). However, the profile com- 



puted with the FALX model in first panel of Fig- 
ure H] for the even isotopes is in good agreement 
with the one given in their paper. 

The profiles in the second panel of Figure HI 
which represent the odd isotopes, also exhibit a 
similar large sensitivity to the choice of model at- 
mosphere. Therefore the combined even-odd iso- 
topes profiles in the third panel are also very sen- 
sitive to the temperature structure. 

For the sake of clarity, let us point out that 
the combined Q/I profiles in the bottom panel of 
Figure [4] differ profoundly from what one would 
obtain from a linear superposition of the corre- 
sponding profiles for the even and odd isotopes 
individually in the two other panels, in propor- 
tion to their isotope ratios. The reason is that the 
combination is highly non-linear, since the lines 
are formed in an optically thick medium (namely 
the radiative transfer effects). While the opacities 
and redistribution matrices are combined in a lin- 
ear way as described by Equations ([2]) and (O, the 
even and odd isotopes blend with each other in the 
radiative transfer process, which makes the com- 
bination as it appears in the emergent spectrum 
highly non-linear. 

The drastic depolarization of Q/I in the line 
core has its origin in the polarizability factor W2 
of the odd isotopes. It is well known that the 
trough like suppression of W2 in the line core for 
Ba II is formed du e to hyperfine structure for the 
odd isotopes (see Istenflol Il997t ). In our radia- 
tive transfer calculations we have a superposition 
(in the proportion of the isotope ratios), of the 
trough-like scattering opacity of the odd isotopes, 
with the peak-like scattering opacity of the even 
isotopes. The shape of the Q/I profile depends 
on the details of radiative transfer and PRD (see 
Equations d2J)— (J6j) ) namely on how these two scat- 
tering opacities non-linearly blend to produce the 
net result for the emergent radiation in the opti- 
cally thick cases. 



5.1.3. Need for the FALX model 

As seen from the last panel of Figure fj] the 
central peak is not well reproduced by any of the 
standard model atmospheres. All the models pro- 
duce a dip at line center. Such a central dip is 
commonly due to the effects of PRD, caused by 
the properties of the type-II frequency redistribu- 
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tion. In the case of Ba n D2, the contribution to 
this central dip comes mainly from the even iso- 
topes, as shown in Figure O The three rows in 
this figure represent the even, odd and combined 
even-odd cases, respectively, for the FALX model 
atmosphere. The first column shows the profiles 
computed with only type-II redistribution, the sec- 
ond column those computed with CRD only. The 
CRD profiles are obtained by setting the branch- 
ing ratios A — and — (1 — e). None of the 
CRD profiles shows a central dip. 

The occurrence and nature o f this central dip 
has b een explored in detail in Holzreuter et al.l 
(|2005l ) , who showed that its magnitude is strongly 
dependent on the choice of atmospheric parame- 
ters. This behavior is also evident from Figure [4] 
The cooler the atmosphere, the smaller is the cen- 
tral dip. The dip is often smoothed out by instru- 
mental and macro-turbulent broadening. How- 
ever, the profiles in Figure 0] have already been 
smeared with a Gaussian function having a full 
width at half maximum (FWHM) of 70 mA. The 
dip could be suppressed by additional smearing, 
but such large smearing would also suppress the 
observed side peaks of the odd isotopes and would 
make the intensity profile inconsistent with the ob- 
served one. The value 70 mA has been chosen to 
optimize the fit, but it is also consistent with what 
we expect based on the observing parameters and 
turbulence in the chromosphere. 

The failure of all the tried standard model 
atmospheres therefore leads us to introduce a 
new model with a modified temperature structure, 
which is cooler than the standard FALX model. 
The details of the new cooler FALX model has 
been given in Section 14.21 This new model at- 
mosphere succeeds in giving a good fit to both the 
intensity and the polarization profiles, as shown in 
Figure [BJ To simulate the effects of spectrograph 
stray light on the intensity and polarization pro- 
files we have applied a spectrally flat unpolarized 
background of 4% of the continuum intensity level 
to the theoretical (I, Q/I) profiles. For a good 
Q/I line center fit, we find that it is necessary to 
include Hanle depolarization from a non-zero mag- 
netic field. Our theoretical profiles are based on 
a micro-turbulent magnetic field of strength -Bturb 
with an isotropic angular distribution. Our best fit 
to the Q/I profile corresponds to a field strength 

Of Bturb = 2 G. 



5.2. The importance of PRD 

The importance of PRD in modeling the 
Ban D2 line has alre ady been demonstrated in 
Faurobert et al. (2009), although by only consid- 
ering the even isotopes. Figure [5] demonstrates 
the importance of PRD for both the odd and the 
even isotopes. As seen from the second column of 
this figure, the Q/I profiles for the odd isotopes 
when computed exclusively in CRD do not pro- 
duce any side peaks, while the profiles computed 
with type-II redistribution exhibits such peaks. 
Also, by comparing the Q/I profiles for the even 
isotopes in the first row, we see that CRD fails 
to generate the needed line wing polarization. A 
comparison between the observed profiles and the 
theoretical profiles based on CRD alone (dotted 
line) and on full PRD (dashed line) for the FALX 
model is shown in Figure [7J While the inten- 
sity profile can be fitted well using either PRD or 
CRD, the polarization profile cannot be fitted at 
all with CRD alone. PRD is therefore essential to 
model the Q/I profiles of the Ba 11 D2 line. 

6. Conclusions 

In the present paper we have for the first time 
tried to model the polarization profiles of the Ba 11 
D2 line by taking full account of PRD, radiative 
transfer, and HFS effects. We use the theory of in- 
state interference developed in PI in combination 
with different atom models representing different 
isotopes of Ba 11 and various choices of model at- 
mospheres. Applications of the well known stan- 
dard model atmospheres FALF, FALC, FALA, 
and FALX fail to reproduce the central peak, and 
instead produce a central dip mainly due to PRD 
effects. We have shown that in the case of Ba 11 
D2 the central dip is reduced by lowering the tem- 
perature of the atmospheric model. We can there- 
fore achieve a good fit to the observed polarization 
profile by slightly reducing the temperature of the 
FALX model. 

In modeling the Ba 11 D2 line we account for the 
depolarizing effects of elastic collisions with hy- 
drogen atoms but neglect the alignment transfer 
between the 2 i-3/2 level and the metastable 2 -Pr/? 



level. It has been shown by iDerouichl ([20081 ) that 
this alignment transfer affects the line center po- 
larization and is needed for magnetic field diagnos- 
tics. The purpose of the present paper is however 
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not to determine magnetic fields but to clarify the 
physics of line formation. 

We demonstrate that PRD is essential to repro- 
duce the triple peak structure and the line wing 
polarization of the Ba n D2 line, but find that the 
line center polarization is very sensitive to the tem- 
perature structure of the atmosphere, which con - 
tradicts the conclusions of lFaurobert et aL (2009), 
who find that the barium line is temperature in- 
sensitive and therefore suitable for Hanle diagnos- 
tics. This contradiction illustrates that a full PRD 
treatment as done in the present paper, including 
the contributions from both the even and odd iso- 
topes, is necessary to bring out the correct tem- 
perature dependence of the line. The large tem- 
perature sensitivity of the Ba II D2 line makes 
it rather unsuited for magnetic-field diagnostics, 
since there is no known straightforward way to 
separate the temperature and magnetic-field ef- 
fects for this line. 

We would like to thank Dr. Michele Bianda for 
useful discussions. 
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Fig. 2. — CCD image showing the I and Q/I spec- 
tra of the Ba n D2 line at \i = 0.1. The ob- 
servations were obtained on June 3, 2008, with 
ZIMPOL-2 at the French THEMIS telescope on 
Tcncrifc. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



Table 1: Wavelengths (A) of the hyperfine transi- 
tions for the odd isotopes of Ba 11 
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Fig. 3. — Top panel: The temperature structure of some of the standard model atmospheres. FALX rep- 
resents the model for which the temperature is reduced by about 300 K over a 700 km range around the 
height of formation of the Ban D2 line. Bottom panel: Comparison between the observed center-to- limb 
variation (CLV) of the continuum intensity and the predictions from different model atmospheres including 
FALX for a wide range of wavelengths from the violet to the IR region of the spectrum. For all the /x values 
the dashed and the dash-triple dotted lines are practically indistinguishable as the models FALX and FALX 
produce nearly identical fits. 
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Fig. 4. — Comparison between the observed Q/I profile and the theoretical profiles for some of the standard 
model atmospheres, separately displayed for the even, odd, and combined even-odd cases. The theoretical 
profiles represent the non-magnetic case and have been smeared with a Gaussian having a full width at half 
maximum (FWHM) of 70 mA to account for instrumental and macro-turbulent broadening. 
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Fig. 5. — Theoretical Q/I profiles computed for the non-magnetic FALX model for the even (first row), 
odd (second row), and combined even-odd (third row) isotopes, with only type-II frequency redistribution 
(first column) and only complete frequency redistribution CRD (second column). A prominent central dip 
is present for the type-II redistribution profiles although they have been smeared with a Gaussian having 
FWHM = 70 mA. 
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Fig. 6. — Fit to the observed profile using the FALX model with and without a micro-turbulent magnetic 
field Bturb- The theoretical profiles have been smeared using a Gaussian with FWHM = 70 mA. 
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Fig. 7. — Comparison between the observed Stokes profiles and the profiles computed with CRD (dotted line) 
and PRD (dashed line) for the FALX model. The theoretical profiles have been smeared with a Gaussian 
having FWHM = 70 mA. The strength of the micro-turbulent magnetic field -Bturb has been chosen to be 
2 G for the PRD and 5 G for the CRD profiles. The dashed line in this figure is same as the dashed line in 
Figure |U 
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